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Abstract 

The increased use of veterinary antibiotics in modern agriculture for therapeutic uses and growth promotion has 
raised concern regarding the environmental impacts of antibiotic residues in soil and water. The mobility and 
transport of antibiotics in the environment depends on their sorption behavior, which is typically predicted by 
extrapolating from an experimentally determined soil-water distribution coefficient (KJ. Accurate determination of 
K d values is important in order to better predict the environmental fate of antibiotics. In this paper, we examine 
different analytical approaches in assessing K d of two major classes of veterinary antibiotics (sulfonamides and 
macrolides) and compare the existing literature data with experimental data obtained in our laboratory. While 
environmental parameters such as soil pH and organic matter content are the most significant factors that affect 
the sorption of antibiotics in soil, it is important to consider the concentrations used, the analytical method 
employed, and the transformations that can occur when determining K d values. Application of solid phase extraction 
and liquid chromatography/mass spectrometry can facilitate accurate determination of rCj at environmentally relevant 
concentrations. Because the bioavailability of antibiotics in soil depends on their sorption behavior, it is important to 
examine current practices in assessing their mobility in soil. 

Keywords: Veterinary pharmaceuticals, Antibiotics, Antimicrobials, Sulfamethazine, Tylosin, Partition coefficient, 
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Introduction 

Veterinary pharmaceuticals (VPs) are physiologically ac- 
tive compounds that are used to protect animals against 
parasites, prevent bacterial infections, and growth promo- 
tion [1-6]. Antibiotics, their synthetic analogues, and syn- 
thetically produced antimicrobials are introduced in both 
therapeutic and medical dosages to the animals through 
medicated feeds, injections, and external application [1]. 
Most animals do not absorb these compounds completely; 
hence antibiotics are excreted in urine and feces as parent 
compounds, and in mixtures with their conjugated metab- 
olites and oxidation/hydrolysis products [1,3-8]. 

The presence of antibiotic residues in animal manure 
that are land-applied to agricultural lands can contamin- 
ate water and soil [9-11]. Consequently, changes in the 
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soil microbial population can occur; the microorganisms' 
ability to degrade contaminants and their role in chem- 
ical cycles, such as nitrification, may be affected signifi- 
cantly [5,7,12-14]. Persistent antibiotics can accumulate 
in the top layers of soil, may leach to the groundwater, 
or can be transported to surface waters [1,15,16]. Sul- 
fonamide antibiotics, the first broad spectrum antibac- 
terial medications [17], are relatively persistent in the 
environment and do not sorb strongly to soil. Thus, sul- 
fonamides have been detected in surface water, ground 
water, soil pore water [18-20], and drinking water [21] 
that have been impacted by agricultural and human 
activities. 

A common parameter used to predict the transport be- 
havior of organic contaminants in soil is the soil-water 
partition coefficient, Kd- The K d values can be directly de- 
termined experimentally, or derived indirectly from the 
octanol-water partition coefficients (K ow ) or by computa- 
tional modeling using free energy calculations. Because 
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the sorption properties of antibiotics affect their mobility 
and ecotoxicology, it is important to recognize that differ- 
ent Kd measurements may provide varying results that 
could potentially lead to large errors in environmental 
models that are used in risk assessment. 

The fate and transport of antibiotics in the environ- 
ment depend on the underlying physical properties of 
the compound such as water solubility, lipophilicity, 
volatility, and sorption potential. Soil can act as potential 
sink, and thus sorption of antibiotics in the solid phase 
can reduce their mobility, reactivity, and bioavailability 
for microbial degradation [22]. In addition, soil proper- 
ties such as organic carbon content, ionic strength, clay 
content, texture, and pH can alter sorption mechanisms 
involved, and the extent of sorption of antibiotics [23]. 
The assumption that sorption occurs solely through 
hydrophobic partitioning to soil organic matter (OM) is 
inappropriate for antibiotics with ionizable groups, when 
electrostatic interactions and hydrogen bonding become 
significant [1]. 

The purpose of this paper is to provide a review of the 
different ways that Kd values are measured, and demon- 
strate how the analytical differences may affect the pre- 
diction of the fate and transport of antibiotics in the 
environment. Specifically, this review will focus on two 
of the most used classes of antibiotics: sulfonamides and 
macrolides. Within these classes, sulfamethazine, tylosin, 
and erythromycin are examined due to their wide use in 
animal related practices, and their variable sorption 
properties [1,24]. 

Review 

Octanol-water partition coefficient (K ow ) and octanol- 
water distribution ratio (D ow ) 

Antibiotic mobility in soil has been traditionally esti- 
mated using the octanol-water partition coefficient 
(K ow ): 

K = \Solute\ octanol 
[Solute] water 

However, K ow only reflects hydrophobic interactions 
and does not accurately account for electrostatic interac- 
tions, surface complexation, hydrogen bonding, cation- 
exchange, or bridging that may vary significantly with 
changes in pH, OM, and ionic strength [1]. While the 
use of K ow in predicting soil sorption behavior of nonpo- 
lar compounds works fairly well, the application of K ow 
for polar or ionizable compounds, such as many antibi- 
otics, may be inaccurate. Additionally, the variety of en- 
vironmental factors (such as soil properties) that can 
affect sorption will complicate the modeling efforts to 
predict sorption and mobility of antibiotics. For ex- 
ample, OM may block interlayer sites of clay minerals 



[25], but such phenomenon is not accounted for by 
K ow values. Thus simply using K ow will result in incor- 
rect assessment of antibiotic fate and transport in the 
environment. 

The pH dependent octanol-water distribution ratio 
D OWj can be used to avoid variations in K ow values result- 
ing from changes in pH. The D ow value considers hydro- 
phobicity and ionogenicity, and is a combination of the 
Kow (of the neutral compound) and the pK a , in which 
the transfer of both neutral and ionized species between 
the aqueous and immiscible phase is accounted for [26]. 
The D ow value does not consider hydrophobicity as the 
sole governing factor that dictates partitioning of neutral 
compounds, but also accounts for the transfer of ion pairs 
and free ions from aqueous to the organic layer [27,28]: 

[nonionized + ionized species] octanol 
[nonionized + ionized species] water 

A relationship between log K ow and log D ow can be de- 
rived for both acidic and basic compounds [29]. For sul- 
famethazine, K ow values between 1.042 and 3.750 are 
reported, while D ow values between 0.427 and 1.950 are 
reported (determined at pH values of 4-8). These K ow 
and D ow values were calculated using Advanced Chemistry 
Development (ACD/Labs) Software VI 1.02 (© 1994-2012 
ACD/Labs). Notably, these values fall in the lower part of 
the wide range of K d values (0.23-30 L/kg) obtained experi- 
mentally, as reported from literature (Table 1). 

The K d partition coefficient 

The soil-water partition coefficient (Kd) is used to de- 
scribe the sorption potential of pollutants and the extent 
that they will move into the ground or surface waters. 
Using K d instead of K ow demonstrates sorption behavior 
with respect to the soil media of interest, and data ex- 
trapolation from the octanol to soil matrices is eliminated. 
The K d value is the ratio between the concentration of the 
compound in soil (C s ) [total concentration, including 
sorbed transformation products] to the concentration of 
the dissolved compound in water (Cw) [1]: 



In the experimental determination of Kd values, it is 
important to accurately measure the concentrations of 
the compounds at environmentally relevant levels for 
both the water and soil components to assure mass bal- 
ance. Kd is typically determined one of two ways: (1) col- 
umn displacement studies in which determination 
occurs from a breakthrough curve at a single location, 
or (2) batch sorption experiments in which multiple 
concentrations are used to construct isotherms by plot- 
ting C s versus C w - 
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Table 1 Sorption coefficients of sulfamethazine 


Soil type 


Kd (L/kg) 


Concentration range used 


Measurement technique 


Reference 


Sandy loam 


0.27-0.77 


1.5, 3.5, 7.5, 10, 15 ppm 


HPLC-UV (X = 275 nm) 


[34] 




0.23-1.22 


0.3 - 20 ppm 


HPLC-UV (X = 254 nm) 


[33] 




9.8-22 


0, 0.25, 0.50, 1 .0, 2.5, 4.0 ppm 


HPLC-UV (X = 275 nm) 


[32] 




0.95-19.53 


17.7, 35.4, 53.1, 70.8, 88.5 ppb 


Liquid scintillation counting 


ED 








(LSC) with 14 C-SMZ 


Method 1 




1.0-7.52 


1, 3, 10, 20, 30, 50, 100, 300 ppb 


LSC with 14 C-SMZ 


ED 










Method 1 


Clay loam 


2.88 


1.5, 3.5, 7.5, 10, 15 ppm 


HPLC-UV (X = 275 nm) 


[34] 




16.55 ± 1.41 


0.012, 0.122, 1.219 ppm 


LSC with 14 C-SMZ 


[31] 


Loam 


1.05-3.91 


0.3 - 20 ppm 


HPLC-UV (X = 254 nm) 


[33] 




3.1-17 


0, 025, 0.50, 1 .0, 2.5, 4.0 ppm 


HPLC-UV (X = 275 nm) 


[32] 




2.83-22.28 


17.7, 354, 53.1, 70.8, 88.5 ppb 


Liquid scintillation counting 


ED 








(LSC) with 14 C-SMZ 


Method 1 




0.9-18.2 


1, 3, 10, 20, 30, 


LSC with 14 C-SMZ 


ED 






50, 100, 300 ppb 




Method 1 




17.10 ± 1.66 


0.012, 0.122, 1.219 ppm 


LSC with 14 C-SMZ 


[31] 


Silty clay loam 


18.58 ± 2.29 


0.012, 0.122, 1.219 ppm 


LSC with 14 C-SMZ 


[31] 


Silt loam 


0.82-2.12 


1.5, 3.5, 7.5, 10, 15 ppm 


HPLC-UV (X = 275 nm) 


[34] 




0.66 - 6.73 


2.5 - 50 uM 


HPLC-UV (X = 254 nm) 


[30] 








LSC with 14 C-SMZ 






206.18 ± 12.09 


0.012, 0.122, 1.219 ppm 


LSC with 14 C-SMZ 


[31] 


Loamy sand 


2.3-30 


0, 0.25, 0.50, 1 .0, 2.5, 4.0 ppm 


HPLC-UV (X = 275 nm) 


[32] 


Sand 


7.52 ± 0.09 


0.012, 0.122, 1219 ppm 


LSC with 14 C-SMZ 


[31] 



Note: Results are soil dependent. ED: experimentally determined. Methodology can be found in the additional files. K ow of sulfamethazine ranges from 
1.042-3.750. K ow values calculated using Advanced Chemistry Development (ACD/Labs) Software V1 1.02 (© 1994-2012 ACD/Labs). 



Experimentally determined K d values reported in 
the literature for a particular compound are highly 
variable even for the same soil type and environmen- 
tal conditions. For example, Tables 1 and 2 list Kj 
values for sulfamethazine, a sulfonamide [30-34] and 
tylosin, a macrolide [2,6,25,35-38], respectively. Cor- 
responding plots were drawn in Figure 1 to clearly 
demonstrate the wide range of their Kd values re- 
ported. For the same type of soil, K d values appear to 
vary widely depending on the concentration ranges 
used to determine I<d- 

In general, Kd values have been obtained using high 
concentrations (in the parts per million range) of antibi- 
otics that are not environmentally relevant. High con- 
centrations are used to enable detection of the desorbed 
portion without sample pre-concentration. However, 
the use of high concentrations of antibiotics for sorp- 
tion experiments can result in anomalies when the Kd 
value is concentration-dependent and exhibits non- 
linearity. For example, the K d value for sulfachloropyri- 
dazine (pk a = 1.88, 5.90) was determined in sandy loam 
(pH 6.0-7.5 and 6.6, respectively) at 1-10 ppb to be 0.9 L/kg, 



while when determined at 1.5 ppm the I<d value was 
8.1 L/kg. These discrepancies in the Kd values pose differ- 
ences in predicting the fate of sulfachloropyridazine; it im- 
plies that at lower concentrations, sulfachloropyridazine is 
considered to have high mobility under the pesticide mo- 
bility classification [6,19], while at higher concentrations 
sulfachloropyridazine has low mobility [6,39] . 

The freundlich sorption constant, K f : 

Because sorption coefficients are not always the same at 
all aqueous concentrations, linear plots are not always 
observed. In the case of tylosin, non-linearity has been 
previously reported [40,41], and therefore all original 
data will be presented as both K d and Kf values. The 
Freundlich constant (sorption coefficient) Kf provides a 
better estimate of partitioning: 



where n, the Freundlich exponent, is a measure of the iso- 
therm nonlinearity. A plot of log C s vs. log C w gives a 
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Table 2 Sorption coefficients of tylosin 



Soil Type 


Kd (L/kg) 


Concentration range used 


Measurement technique 


Reference 


Sandy Loam 


58.1 - 148.0 


500 ppb 


LC-MS 


[37] 








LC-fluorescence 






92 


1000 ppm 


ELISA 


[35] 




101.02-13961.00 


10, 100, 200, 1000 ppb 


HPLC-MS 


ED 








(ion trap) 


Method 2 




6737-33871 


1, 5, 10, 100, 200, 500, 1000 ppb 


HPLC-MS/MS 


ED 










Method 3 




42 


5 ppm 


HPLC-UV 


[36] 








(X = 290 nm) 




Clay Loam 


66 


1000 ppm 


ELISA 


[35] 




156 ± 8 


3 ppm 


HPLC-MS (SIM) 


[6] 




1.76 - 6.19 


0.5, 5, 50 ppm 


HPLC-UV 


[25] 








(X = 290 nm) 






65 


5 ppm 


HPLC-UV 


[36] 








(X = 290 nm) 




Loam 


59.35-1 1 76.00 


10, 100, 200, 1000 ppb 


HPLC-MS (ion trap) 


ED 










Method 2 




1684-172480 


1, 5, 10, 100, 200, 500, 1000 ppb 


HPLC-MS/MS 


ED 










Method 3 




5.77 - 12.4 


0.5, 5, 50 ppm 


HPLC-UV 


[25] 








(X = 290 nm) 




Loamy Sand 


8.3 ± 1 .2 


500 ppb 


LC-MS 


[37] 








LC-fluorescence 






8.9 ± 0.4 


3 ppm 


HPLC-MS (SIM) 


[6] 


Sand Soil 


10.8 ± 0.7 


500 ppb 


LC-MS 


[37] 








LC-fluorescence 






24 


5 ppm 


HPLC-UV 


[36] 








(X = 290 nm) 




Clay%: 0.2 - 51.6% 


1 0.4-387.0 


3 - 7.5 ppm 


HPLC-MS (SIM) 


[38] 




(Average = 1 29.5) 








Clay%: < 3 - 69% (Kj increases with clay %) 


2.23-5520 


5 - 43 umol/L 


HPLC-UV, 


[2] 




(Tylosin A) 




(X = 280 nm) 






547-4745 


5 - 43 umol/L 


HPLC-UV, 


[2] 




(Tylosin D) 




(X = 280 nm) 






597-6520 


5 - 43 umol/L 


HPLC-UV, 


[2] 




(Tylosin A-Aldol) 




(X = 280 nm) 





Note: Results are soil dependent. ED: experimentally determined. Methodology can be found in the additional files. K ow ranges for tylosin A and it's metabolites: 
tylosin A: 0.552-32.6587; tylosin B: 1.535-78.343; tylosin C: 0.676-41. 495;tylosin D: 0.309-17.947. K ow values calculated using Advanced Chemistry Development 
(ACD/Labs) Software V1 1.02 (© 1994-2012 ACD/Labs). 



linear isotherm with a slope equal to n and a y-intercept 
equal to log Kf. If the value of Kf approaches the value of 
Kj, the Freundlich exponent, n, is equal to 1, and sorption 
is linear. If n is greater than 1, the sorption coefficient 
increases as the amount of compound sorbed on the 
solid phase increases; this indicates that the presence 



of sorbed compounds on the solid induces further 
sorption of additional compounds. If n is less than 1, 
sorption coefficient decreases when the amount of 
compound sorbed is increased; this indicates that the 
presence of sorbed compounds hinders further sorp- 
tion [29]. 
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Figure 1 Box plots of Kj values for sulfonamides and macrolides reported in literature. The sulfonamides (left) include sulfamethazine and 
sulfachloropyridazine, and the macrolides (right) includes tylosin and erythromycin. pH values range from 5.2-7.5 when reported. Soil types 
include loamy and sandy loam, clay loam, loam, loamy sand, and silt loam. The high variation of K d values found in literature is illustrated here. 
The upper and lower boundaries of the box represent the 75 th and 25 th percentile respectively. The middle line indicates the median value, and 
the whiskers indicate the maximum and minimum values. 



An alternative to K d : normalizing with organic carbon, K oc 

Experimental determination of Kj values can be cost- 
prohibitive and time-consuming because one has to 
measure Kj at various conditions (e.g. different soil 
types, pH values, and organic and ionic strengths). 
When K d is normalized to the organic carbon content 
of the soil, the organic carbon normalized sorption co- 
efficient K oc is obtained [1]: 



However, mechanisms other than hydrophobic interac- 
tions are not accurately accounted for when normalization 
is performed using organic carbon content [1]. The differ- 
ences between K oc and Kd are observed in literature. 
Rabolle and Spliid [37] reported Kd and K oc values ranging 
from 8.3-128 L/kg and 553-7988 L/kg, respectively, for 
tylosin in 4 different soils. Lertpaitoonpan et al. [33] re- 
ported Kd values for sulfamethazine for 5 different soils at 
varying pH between 0.23-3.91 L/kg, and K oc values be- 
tween 30.4-139.7 L/kg. In both cases, the antibiotics have 
higher K oc values, which would suggest that the com- 
pounds are less mobile than their Kd values would indi- 
cate. Thus, while normalizing partition coefficients may 
help reduce variation between samples, it cannot be uni- 
versally applied to all antimicrobials, particularly those 
that have ionizable functional groups. 



Case studies: sorption behavior of sulfonamides and 
macrolides in sediment 

Macrolides and sulfonamides are commonly used anti- 
biotic classes in livestock. Approximately 165800 kg of 
tylosin (a macrolide), 18660 kg of sulfamethazine, and 
19400 kg of sulfathiazole are used annually in the United 
States for growth promotion, prevention, and therapy 
[42]. Our laboratory conducted sorption experiments for 
sulfamethazine and tylosin under varying pH, OM con- 
tent, and ionic strengths using loam and sandy loam 
sediments. A study by Kim et al. [43] found sulfametha- 
zine, erythromycin-hydrochloride, and tylosin in agricul- 
tural soils at concentrations of 9.1, 30.7, and 19.6 ug/kg, 
respectively. Therefore, sorption tests were performed 
using aqueous concentrations between 1-1000 ug/L 
prior to partitioning in order to mimic environmentally 
relevant concentrations of these antibiotics. Details re- 
garding the methodology used to perform these batch 
experiments can be found in Additional file 1. 

The pH-dependence of antibiotic sorption is critical, be- 
cause many pharmaceuticals have acid-base properties 
resulting in changes in the overall net charge of the mol- 
ecule as ammonia concentration in manure changes [6]. 
These factors can alter the distribution between the aque- 
ous and solid phase, particularly for ionizable compounds 
[23]. Changes in soil pH can also affect surface charge and 
cation exchange capacity of the soil [1]. Ionic strength vari- 
ations can lead to changes in pH, and cause electrostatic 
competition between ions present in the solution and the 
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analyte of interest [6,44-46]. This study makes use of sedi- 
ments that have similar OM content but have different 
fractions of sand, silt and clay. Most of the study conditions 
render a percentage of the compound in its ionized form, 
and due to the dependence on ionic strength, the antibi- 
otics in the cationic form show increased sorption. How- 
ever, sediment buffering capacity must be considered. The 
higher clay content in the loam sediment has a weaker 
buffering capacity relative to OM [47]. The sandy loam can 
more readily adjust the pH closer to the original pH, and 
therefore antibiotic sorption in sandy loam is less affected 
by changes in pH. The water solubility of the antibiotics in- 
creases with increase in dissolved OM content [48], which 
in turn results in increased mobility of antibiotics in soil 
[1]. Thus, it is important to understand how the Kj changes 
for each antibiotic when OM is present in the system. 

Sulfamethazine 

Sulfonamides, or sulfa drugs, are synthetic antimicro- 
bial agents containing the sulfonamide functional group 
(-RS0 2 NH 2 ) [10]. Sulfonamides are mobile antibiotics and 
their speciation changes with pH. A common sulfonamide 
antibiotic is sulfamethazine (pk a : 1.62, 7.91), and its Kd 
values for various soil types reported in literature are pre- 
sented in Table 1, and compared with the Kd values ob- 
tained experimentally from our laboratory. 

Effects of pH on sulfamethazine sorption 

The sorption isotherm we determined for sulfametha- 
zine (Figure 2) illustrates that Kd generally decreases 
with increase in pH for both loamy sand and loam sedi- 
ments. This sorption behavior is consistent with the 
changes in the fraction of ionization of sulfamethazine 
as it converts from its cationic form to the neutral and 
anionic forms (See Figure 3). Positively charged species 
are electrostatically attracted to the negatively charged 
soil surface, and therefore a higher Kd is observed at pH 
below 5 (Table 3) [2,30]. Despite the presence of a small 
fraction of negatively charged sulfamethazine at pH 7, 
cation bridging does not appear to play a significant role 
in the sorption of sulfamethazine because sulfonamides 
interact primarily with soil organic matter via hydropho- 
bic interactions [49]. This behavior of sulfonamides is in 
contrast with tetracycline and fluoroquinolone family of 
antibiotics that interact with soils primarily through cat- 
ion exchange, surface complexation and cation bridging 
sorption mechanisms. 

Sulfamethazine sorption (Table 3) trends towards lin- 
ear isotherms in the sandy loam (n pH s = 0.916, n pH 7 = 
0.853, and n pH 9 = 1.01) at the three pH values tested. 
Sorption in loam exhibits some non-linearity (n pH4 = 
0.885, n pH 6.9 = 0.822, and n pH 8.2 = 0.708). The n values 
are less than 1 which signifies that the sorption coeffi- 
cient decreases when the amount of compound sorbed 



is increased, indicating that the presence of sorbed com- 
pounds hinders further sorption of antibiotic [29,50]. 

Effects of ionic strength on sulfamethazine sorption 

Sulfamethazine showed a slight decrease in sorption 
when ionic strength was increased from 50 mM to 250 
mM (Table 4). The small decrease in sulfamethazine 
sorption may be attributed to a slight change in pH 
brought about by increase in ionic strength, and a pos- 
sible change in interfacial potential between the negative 
sediment surface and the partially charged sulfametha- 
zine [44-46]. The negatively charged soil surface reduces 
the sorption of anionic organic compounds [6,51]. 

Effects of organic matter on sulfamethazine sorption 

Fan et al. [31] found that the sorption correlation of 
sulfamethazine with OM is confounded by soil pH. 
Thiele-Bruhn and Aust [52] observed that when elec- 
trostatic competition were eliminated through use of 
an acidic pig slurry matrix, the sorption of sulfon- 
amides decreased. Decreased antibiotic sorption may 
be attributed to association of sediments with OM 
from manure components (ammonia-N-containing sol- 
uble hydrocarbons such as amino acids urea [53,54], and 
N-heterocyclic hydrocarbons such as pyrrols, methylin- 
dols, and nitrogen bases [55]). The interaction of OM 
with soil can block access of antibiotics to interlayer 
sorption sites in soil [25,52,56]. In our study, we found 
no consistent trend with changes in humic acid (HA) 
concentrations (Table 5). These tests may have been 
complicated by the presence of both dissolved and sus- 
pended HA in solution. Suspended HA provides sites 
where additional partitioning can occur. Increased amounts 
of dissolved OM can cause antibiotics to desorb from soil, 
and increased association of antibiotics with dissolved OM 
can facilitate transport in the environment [1,16,57]. Fur- 
thermore, any anionic sulfamethazine may be repelled by 
the increased surface charge occurring from the dissolved 
OM [30]. 

Tylosin 

Macrolides, which are mainly active to Gram-positive 
bacteria, inhibit ribosomal protein synthesis. Their activ- 
ity stems from the presence of the macrolide ring, a 
large lactone ring to which one or more deoxy sugars 
are attached [58]. A case study on the soil sorption of 
tylosin antibiotic, which belongs to the macrolide class, 
is presented below. 

Effects of pH on tylosin sorption 

Tylosin sorption (pK a : 7.20, 12.44, 12.93, 13.36, 13.94, 
and 15.01; assignments of pKa values in the molecule 
are shown in Scheme 1) strongly depends on the pH, as 
well as on the surface area, clay content, and cation- 
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Figure 2 Sulfamethazine sorption isotherms. Top: Sulfamethazine sorption isotherms in sediment at low, neutral, and high aqueous pH. Left: 
sandy-loam and Right: loam Bottom: Tylosin sorption isotherms in sediment at low, neutral, and high aqueous pH. Left: sandy-loam and 
Right: loam 



exchange capacity of the soil [2]. Since tylosin is water 
soluble (5 mg/mL) and has high molecular weight, it is 
unlikely that sorption occurs through penetration of 
soil micro pores [25]. Several studies have reported that 
the K d values for tylosin increase with decreasing pH 
[6,25,38,59]. The same pH effects on the sorption behav- 
ior of tylosin were observed in the studies conducted in 



our laboratory, as shown in Figure 2 and in Table 3. 
Tylosin sorption increased in both loam and sandy loam 
sediments when the pH of the sediment-aqueous system 
was decreased. The increased sorption of tylosin at 
pH 5, relative to its sorption at pH 7 and 9 can be ex- 
pected due to the shift in tylosin speciation towards 
the positively charged species, resulting in increased 




Figure 3 Distribution of sulfamethazine species by pH. Sulfamethazine is predominantly cationic below pH 1.62, neutral between pH 
1.62 - 7.91, and anionic above pH 7.91. The chemical structures that represent the highest fraction of species is shown above the curve. 
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Table 3 Sulfamethazine and tylosin partitioning with changes in pH 



Sandy loam Loam 



Antibiotic 


PH 


K d (L/kg) 


R 2 


K f 


R 2 


Kd (L/kg) 


R 2 


K f 


R 2 


Sulfamethazine 


5 


6.9 ± 0.5 


0.996 


8.49 


0.997 


18± 1 


0.998 


24.3 


0.999 




/ 


5.1 ±0.3 


0.997 


7.54 


0.996 


5.3 ±0.7 


0.985 


10.3 


0.998 




9 


0.9 ±0.3 


0.954 


0.994 


0.980 


0.9 ±0.1 


0.990 


1.98 


0.995 


Tylosin 


5 


3x1 0 4 ± 1x10 4 


0.933 


6.28x1 0 4 


0.967 


1.7x10 s ±4x10 4 


0.952 


2.02x10 5 


0.933 




7 


1.5x10 4 ±6x10 3 


0.899 


1.43x10 4 


0.950 


2.2x1 0 4 ± 7x1 0 3 


0.933 


1.66x10 4 


0.867 




9 


7x10 3 ±3x10 3 


0.872 


2.21x10 3 


0.962 


1.7x10 3 ±9x10 2 


0.808 


1.26x1 0 3 


0.938 



Note: Values were experimentally determined. Errors represent standard deviation. Methodology can be found in the additional files. Sulfamethazine data was 
performed using method 1 and tylosin data with method 2. Sulfamethazine batch experiments were performed with concentrations of 1, 3, 5, 10, 20, 30, 50, 100, 
and 300 ng/mL, and tylosin batch experiments were performed with concentrations of 1, 5, 10, 100, 200, 500, and ng/mL. 



electrostatic attractions to the negatively charged sedi- 
ment surface [25]. 

Our experimental values indicate that tylosin sorption 
(Table 3) is linear over 3 orders of magnitude in the sandy 
loam sediment at pH 5 and 7 (n pH s = 0.993, n pH 7 = 1.05). 
However, at pH 9, n pH 9 = 1.22; this greater than unity 
value indicates a non-linear sorption behavior that can be 
attributed to the presence of sorbate molecules inducing 
further sorption [29]. In other words, the presence of the 
sorbed tylosin results in further sorption of the antibiotic 
in soil. Similarly, the loam sediment exhibits some non- 
linearity (n pH 5 = l-15, n pH 7=1.18, n pH 9 = l-19), but to a 
lower extent. Thus, electrostatic forces dominate the sorp- 
tion model. 

Effects of ionic strength on tylosin sorption 

Literature suggests that tylosin sorption decreases with 
increase in ionic strength due to the consequent change 
in pH, and as a result of competition between the elec- 
trolyte cations and the positively charged tylosin species 
for negatively charged sorbent [6] . However, ionic strength 
experiments performed in our laboratory at a constant pH 
of 7, with tylosin in mostly neutral form, showed a reverse 
trend (Table 4). Instead, at pH 7, tylosin sorption in- 
creased with increase in ionic strength. This sorption be- 
havior may be attributed to the presence of hydrated 



cations in the solution (Ca 2+ , Na + ) that may act as proton- 
donors, which can protonate the tertiary amine in the 
tylosin molecule and enhance its sorption properties at 
higher ionic strengths. Yong-Hak et al. [60], observed that 
the tertiary amine group of erythromycin can become pro- 
tonated, and that clay surfaces can facilitate this with their 
proton supplying power. Alternatively, hydrated cations 
that adsorb on the negatively charged soil can provide 
hydrogen bonding as an important sorption mechanism 
for tylosin because of several OH groups present in the 
molecule. 

Effect of organic matter on tylosin sorption 

It was expected that the Kd values for tylosin would de- 
crease in the presence of OM due to increased solubility. 
However, studies in our laboratory demonstrated higher 
Kd values with increased OM (represented as humic 
acids) using 10 ppm tylosin in sandy loam, and un- 
changed values in loam (Table 5). Similar to sulfametha- 
zine, complications could arise from the presence of 
both dissolved and suspended OM within the solution. 
This complexity can be observed in the change in K f 
values with increasing tylosin concentration (Table 5). 
Likewise, differences in the sorbates can also influence 
sorption, as was observed in the sorption of tylosin to 
the two sediment types used in our laboratory study. 



Table 4 Sulfamethazine and tylosin partitioning with changes in ionic strength 



Sandy loam 



Loam 



Antibiotic 


Ionic strength (mM) 


Kd (L/kg) 


R 2 


Kf 


R 2 


K d (L/kg) 


R 2 


Kf 


R 2 


Sulfamethazine 


30 


10.3 ±0.5 


0.999 


26 


0.992 


5.5 ± 0.8 


0.991 


46 


0.978 




50 


10± 1 


0.994 


4/ 


0.991 


7.5 ± 0.5 


0.997 


59 


0.960 




250 


7±3 


0.913 


54 


0.974 


2.8 ± 0.3 


0.995 


84 


0.917 


Tylosin 


30 


500 ±45 


0.969 


1640 


0.969 


1 90 ± 30 


0.984 


2689 


0.919 




50 


2.9x1 0 3 ± 1x1 0 2 


0.999 


4411 


0.980 


600 ± 300 


0.886 


3109 


0.972 




250 


1.4x10 4 ±7x10 3 


0.927 


17200 


0.935 


9.6x1 0 3 ± 7x1 0 2 


0.998 


5555 


0.937 



Note: Values were experimentally determined. Errors represent standard deviation. Methodology can be found in the additional files. Sulfamethazine data v 
performed using method 1 and tylosin data with method 2. Sulfamethazine batch experiments were performed with concentrations of 17.7, 35.4, 53.1, 70.£ 
ng/mL, and tylosin batch experiments were performed with concentrations of 10, 100, 500, and 1000 ng/mL 
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Table 5 Sulfamethazine and tylosin partitioning with changes in organic strength 



Sandy loam Loam 



Antibiotic 


Humic acid 
(mg/L) 


Kd 
(L/kg) 


R 2 


K f 


R 2 


K d 
(L/kg) 


R 2 


K f 


R 2 


Sulfamethazine 


1 


9 ± 2 


0.974 


31 


0.910 


5 ± 2 


0.918 


70 


0.913 




10 


3 ± 2 


0.727 


67 


0.947 


4 ± 2 


0.878 


84 


0.827 




50 


12 + 3 


0.960 


50 


0.910 


6 ± 2 


0.961 


86 


0.958 


Tylosin 


1 


60 ± 30 


0.869 


1.43x10 s 


0.873 


50 ± 10 


0.960 


1.66x10 s 


0.983 




10 


90 ± 20 


0.976 


4.41 x10 3 


0.940 


50 ± 50 


0.715 


1.77x10 4 


0.829 




50 


140 ± 70 


0.936 


8.1 8x1 0 4 


0.996 


40 ± 30 


0.714 


1.19x10 s 


0.975 



Note: Vaiues were experimentally determined. Errors represent standard deviation. Methodology can be found in the additional files. Sulfamethazine data was 
performed using method 1 and tylosin data with method 2. Sulfamethazine batch experiments were performed with concentrations of 17.7, 35.4, 53.1, 70.8, 88.5 
ng/mL, and tylosin batch experiments were performed with concentrations of 10, 15, 30, and 50 mg/L. 



Sorption of tylosin metabolites 

Tylosin and other antibiotics may interconvert between 
multiple chemical forms depending on environmental 
conditions as shown in Figure 4. Tylosin A and its related 
compounds are stable from pH 4-9 [61]. Metabolism of 
tylosin by livestock results in the excretion of tylosin A, B, 
D, and dihydrodesmycosin metabolites [2,62,63]. As the me- 
tabolites retain different degrees of bioactivity (TA = 100%, 
relative, TB = 83%, TD = 35%, dihydrodesmycosin = 31%) 
[2], it is important to consider the speciation of tylosin 
present in the environment. Tylosin A, D, and tylosin 
A-Aldol have been found to exhibit similar sorption charac- 
teristics [2]. However differences in sorption behavior be- 
tween tylosin A (K ow : 0.552-32.659) and its hydrolysis 

pK a =12.93 



OH 




pK a =7.20 

Scheme 1 The macrolide, tylosin. Chemical structure and pK a 
values are shown. 



product tylosin B (spiramyycin, K ow : 1.535-78.343) may be 
significant. Tylosin B results from the hydrolysis of tylosin 
A which involves a loss of the mycarose ring attached at 
position 4 of the 16-membered lactone ring. With this loss, 
the hydrophilicity of tylosin increases. This can alter tylo- 
sin's potential to sorb to soil through hydrophobic interac- 
tions, and tylosin B can potentially be more mobile in the 
environment. Therefore, due to the varying properties of 
the different forms of tylosin, including tylosin A, B, C, and 
D it may not be appropriate to use only one Kd value for 
risk assessment of tylosin. Rather, Kd values should be ob- 
tained for all forms possible under the expected conditions. 

To date, studies on the environmental fate of tylosin A 
degradation products are very limited [2,24]. Our labora- 
tory determined sorption differences between tylosin A 
and B in loam sediment at an initial concentration of 
50 ug/mL equilibrated between sediment and aqueous 
phase for 24 h. The amount of tylosin remaining in the 
aqueous phase was determined by liquid chromatography 
coupled to an ion trap mass spectrometer (LC-MS) fol- 
lowing concentration by solid phase extraction (SPE). 
The methodology used follows that in method 2 of the 
additional file 1. It was found that 53% tylosin A and 
39% tylosin B were sorbed in the loam. However, these 
results may be complicated by the hydrolysis of tylosin 
over the equilibration time period and during the SPE 
process. A study by Ali et al. [64] observed a decrease 
in sorption with decreasing pH. This result is contra- 
dictory to what is found in most of the studies reported in 
the literature. The decrease may be associated with the de- 
crease in tylosin A due to the formation of tylosin me- 
tabolites. The conversion of the parent compound to 
metabolites and the interconversions that occur under 
varying conditions are challenges associated with ana- 
lyzing degradation products. 

Conclusions 

It is not possible to determine the fate and mobility of 
antibiotics and antimicrobials in the environment with 



Wegst-Uhrich et at. Chemistry Central Journal 2014, 8:5 
http://journal.chemistrycentral.eom/content/8/1/5 



Page 1 0 of 1 2 




Tylosin D (relomycin) Dihydrodesmycosin 

Figure 4 Degradation products of tylosin. Under environmental conditions, tylosin A can degrade to desmycosin, relomycin, dihydrodesmycosin, 
and tylosin A-Aldol. Tylosin A, relomycin, dihydrodesmycosin.desmycosin, and additional unknown degradates are present in swine excreta [62,63]. 



K ow values alone. The variability in mobility, as demon- 
strated by K d and K f values due to environmental factors 
such as pH, ionic strength, and organic strength as well 
as the multiple chemical functions of the molecule are 
large. In the case studies presented here, sulfamethazine 
was found to be very mobile in sandy loam and loam 
sediments, while tylosin is very immobile in both sedi- 
ments. It is possible that degradates may have a higher 
mobility than the parent compounds. It is also import- 
ant to determine sorption coefficients of antibiotics at 
environmentally-relevant concentrations. To achieve 
this, highly sensitive analytical techniques must be 
used, including the use of radiolabeled compounds. 

A decrease in solution pH resulted in an increase in 
sorption of the cationic forms of antibiotics suggesting 
that electrostatic forces are the favored sorption mech- 
anism of sulfamethazine and tylosin. As with other 
known pharmaceuticals, ionization of these compounds 
at the conditions considered was shown to favor the 
sorption of compounds. A cation-exchange mechanism 
can also be envisioned based on the results of ionic 
strength experiments where ions compete with charged 
species for sites on the soil. Organic matter dependence 
of I<a appears to be concentration-dependent, where low 
antibiotic concentrations result in higher soil sorption, 
and higher antibiotic concentrations result in lower soil 
sorption. It is observed that sorption mechanisms are 
much more complex than simple hydrophobicity and 
hydrogen bonding, and should also consider van der 



Waals and electrostatic interactions, as well as cation- 
exchange, competition and bridging. Additionally, the 
properties of the sorbent also affect the sorption process. 
Differences in clay content alone provide notable changes 
in K d values. Finally, sorption of antibiotics in soil, ma- 
nure, and biosolids can be microbially-mediated, and may 
result in degradation or possibly irreversible binding onto 
manure solids with time [25,40,62,65]. Fate and transport 
studies should take into account not only Kd values for 
the parent compounds, but also those of the transform- 
ation products formed during biotic and abiotic processes 
in soil. 

Additional file 



Additional file 1: Methodology for Sorption Batch Experiments. 
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